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Abstract 

This study investigated the Navier-Stokes 
computations of the surface heat transfer coefficients of a 
transition duct flow. A transition duct from an axisymmetric 
cross section to a nonaxisymmetric cross section, is usually 
used to connect the turbine exit to the nozzle. As the gas 
turbine inlet temperature increases, the transition duct is 
subjected to the high temperature at the gas turbine exit. 
The transition duct flow has combined development of 
hydraulic and thermal entry length. The design of the 
transition duct required accurate surface heat transfer 
coefficients. The Navier-Stokes computational method 
could be used to predict the surface heat transfer coefficients 
of a transition duct flow. The Proteus three- 
dimensional Navier-Stokes numerical computational code 
was used in this study. The code was first studied for the 
computations of the turbulent developing flow properties 
within a circular duct and a square duct. The code was then 
used to compute the turbulent flow properties of a transition 
duct flow. The computational results of the surface pressure, 
the skin friction factor, and the surface heat transfer 
coefficient were described and compared with their values 
obtained from theoretical analyses or experiments. The 
comparison showed that the Navier-Stokes computation 
could predict approximately the surface heat transfer 
coefficients of a transition duct flow. 

Nomenclatures 

C f skin friction factor, t^I/ 2 p^U^) 

D hydraulic diameter 

H heat transfer coefficient, - TJ, Btu/ 

(hr-ft 2 °R) 

h dimensionless heat transfer coefficient, 

HR/K 

K heat conductivity 


k turbulence kinetic energy 

L domain along X-direction 

n normal direction 

Pr Prandtl number, 0.71 

p static pressure 

Q w heat transfer rate 

R circular duct radius 

Re b bulk Reynolds number, p^U^D/p^ 

T temperature 

U,V,W velocity components 

X,Y,Z Cartisian coordinates 

X,r,0 cylindrical coordinates 

Y + dimensionless distance from wall surface, 

n VTwPw A^w 

Px’Py’Pz £rid packing parameters 
E turbulence kinetic energy dissipation rate 

x w shear stress 

Subscripts 

c center-line condition 

d downstream location 

s,w surface/wall condition 
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u upstream inlet location 

oo inlet condition 

o,«> inlet stagnation condition 

Introduction 

Nonaxi symmetric exhaust nozzles are often employed 
on the modem aircraft propulsion systems to improve the 
aircraft performance. A transition duct, from an 
axisymmetric cross section to a nonaxi symmetric cross 
section, is needed to connect the turbine exit to the nozzle. 
In applications, the incoming flow to the transition duct is 
turbulent and subsonic. Existing research 1,2 has 
experimentally explored the aerodynamics of a circular to 
rectangular transition duct flow. Their emphases were on 
the measurements of the total pressure and the shear 
stresses. As the gas turbine inlet temperature increases, the 
transition duct is subjected to the high temperature at the 
gas turbine exit. The transition duct flow has combined 
development of hydraulic and thermal entry length. The 
temperature distribution along the duct surface is an 
important factor in the design of the transition duct. The 
attainment of the duct surface temperature requires the 
accurate predictions of the surface heat transfer coefficients. 

A method of using a liquid crystal-heater composite 
sheet for heat transfer research was developed at the 
NASA Lewis Research Center. 34 The method was 
successful to resolve the heat transfer coefficients of a 
square to rectangular transition duct flow. 5 Recently, this 
method was used to measure the heat transfer coefficients 
of a circular to rectangular transition duct flow. 6 

Numerical computational methods of the Navier- 
Stokes equations and the energy equation could be used to 
compute the velocity and the temperature variations within 
a flow field. Using the advanced computational grid 
generation technique, such as the Gridgen Routines, 7 one 
could optimize the computational grid for complicated 
flow geometries such as the internal flows within transition 
ducts. The turbulence modeling technique has also been 
improved 8 for the Navier-Stokes computations of turbulent 
flows. With proper turbulence model and computational 
grid in the Navier-Stokes computation, it could be used to 
compute the temperature variation within the near wall 
region of a turbulent transition duct flow. Based on the 
near wall temperature variations, we could calculate the 
surface heat transfer coefficients. 

The objective of this study was to investigate the 
capability of the Navier-Stokes code to predict the heat 
transfer coefficients at the wall surface of an experimental 
transition duct. 1 The transition duct, Fig. 1, changed its 
circular cross section at the upstream locations to 
rectangular cross section at the downstream locations. The 


axis of the duct was in line with the upstream flow 
direction. The upstream flow was turbulent and had a 
Mach number of 0.34. The Reynolds number of the 
upstream flow was 1 .6x 1 0 6 . The Proteus three-dimensional 
Navier-Stokes numerical computational code 9 was used 
for the present computational study. The Gridgen software 
program was used to generate the computational grid. The 
unsteady, compressible, and turbulent Reynolds averaged 
Navier-Stokes equations and the energy equation were 
solved numerically for the steady state flow properties of 
the transition duct flow. 

In this report, we describe a preliminary study of 
using the Proteus code to predict the heat transfer 
coefficients of the transition duct flow. The code was first 
studied for its capabilities to calculate the turbulent internal 
flow properties of a circular duct and a square duct. These 
served as baseline verifications of the computational code. 
The code was then used to compute the properties of the 
transition duct flow. The computed heat transfer 
coefficients were compared with their measurements to 
investigate the capability of Navier-Stokes computation 
of the transition duct surface heat transfer coefficients. 

Computation 

The Proteus Navier-Stokes numerical computational 
code was used in this study. The theory, the numerical 
schemes, and the usage of the code could be found in 
Refs. 9 to 1 1 . The code solved the unsteady, compressible 
Reynolds-averaged Navier-Stokes equations and energy 
equation for turbulent flow properties. The Baldwin- 
Lomax eddy viscosity formula 12 and the Chien’s k-e 
turbulence model 13 were coded for turbulence modeling. 
The code solved the governing equations in either the 
Cartesian or the cylindrical coordinate systems. The 
equations were solved by marching in time using the 
generalized time differencing of Beam and Warming. 14 
The second order central differencing was used for all 
spatial derivatives. Nonlinear terms were linearized 
using second-order Taylor series expansions. The finite 
difference equations were solved using an altemating- 
direction-implicit technique. 15 A set of coordinates to 
specify the computational grid was required for a particular 
flow geometry. These coordinates could be created by a 
separate grid generation code and stored in an unformatted 
file for the Proteus code to read. A variety of the boundary 
conditions of mean velocity components, temperature, 
pressure and turbulence for the computations were built 
into the code including (1) specified values and/or 
gradients, and (2) linear extrapolation. Thus, the Proteus 
Navier-Stokes code was suitable for the present heat 
transfer computations of the internal duct flows with 
combined development of hydraulic and thermal entry 
length. 
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In this study, the Proteus code was used to calculate 
respectively the properties of the developing turbulent 
flows within (a) a circular duct, (b) a square duct and (c) a 
circular to rectangular transition duct. The present 
computation used most of the default numerical 
computational schemes in the existing code. The second 
derivative viscous terms were included in the computation 
and it did not use the thin-layer option. The computations 
also included the constant coefficient models of fourth- 
order explicit artificial viscosity and second-order implicit 
artificial viscosity. The coefficients were all 1 .0 for the 
fourth-order viscosity model and all 2.0 for the second- 
order viscosity model. The Euler implicit method and a 
local time step based on the CFL criteria for explicit 
method were utilized in the computations. A CFL value of 
10 was used in the computation with the Baldwin-Lomax 
eddy viscosity formula for turbulence modeling. A small 
CFL value of 1 was used in the computation with the 
Chien’sk-e two-equation model for turbulence modeling. 
The convergence was based on the averaged absolute 
value of the residual for each conservation equation. The 
convergence was assumed when the maximum residual 
was less than 4X10A The boundary conditions and the 
initial conditions were specified for each internal duct 
flow computations. The effect of the computational grid 
and the turbulence model on the computations were first 
studied in the circular duct flow and the square duct flow 
computations. These computations were served as the 
code verifications. The code was then used to compute the 
transition duct flow surface heat transfer coefficients with 
selected grid configurations and turbulence model. The 
grid configurations, the initial conditions, the boundary 
conditions and the turbulence model for each flow 
computation were described in the following section. 

Circular Duct Flow 

The hydrodynamics and thermodynamics of acircular 
duct flow is often used to simulate the internal flow 
properties of noncircular duct. However, existing heat 
transfer studies of a circular duct flow 16 are limited to the 
case with a fully-developed velocity profile condition. In 
practice, the hydrodynamics and thermodynamics flow 
fields developed simultaneously in the inlet region. The 
computation of the heat transfer properties of the 
developing circular duct flow was considered here. Figure 2 
shows a sketch of the flow configuration. The duct has a 
constant circular cross section with a radius of R = 0.36 ft. 
We used the cylindrical coordinate system for the 
computation. 

Boundary Conditions . To simulate the flow field with 
combined hydraulic and thermal entry length, uniform 
mean velocity and mean temperature were assigned to be 
the boundary conditions at X/R = 0. Near the duct inlet 


(X/R = 0), the static pressure was assumed to vary linearly 
along the X-direction. The downstream boundary of the 
computational domain was chosen to be at X/R = 1 00 (L = 
100R). A fully developed flow existed 17 near the 
downstream boundary location. Thus, the X-direction 
gradients of the mean velocity components and the static 
temperature were zero at X/R = 100 location. Uniform 
static pressure was specified at the downstream boundary 
locations. The computations were performed within a 
quadrant of the duct geometry and symmetric boundary 
conditions were assigned at 0 = 0° and 0 = 90° locations. 
The r-direction gradients of the flow properties were zero 
at the duct center line locations. Nonslip conditions were 
the velocity boundary conditions at the duct wall surface. 
The normal pressure gradient was set to zero at the duct 
wall surface. The duct wall surface temperature was 
assumed to be uniform and had a value of 1 . 1 T^. 

The Baldwin-Lomax eddy viscosity formula was 
used for turbulence modeling in the computation. This 
algebraic turbulence model was modified for an internal 
duct flow computations. 9 This turbulence model does not 
require the turbulence boundary conditions for its usage in 
the computation. 

Computational Grids . The grid dimensions were 81 
along the X-direction, 31 along the 0-direction, and 50 
along the r-direction. The grid points at each X location 
were located in the duct cross section. The grid 
configuration was kept the same at every duct cross 
section. The following equation, 

X / L = (p x + 1) - (P x - l){[(p x + l)/(p x - 1)] 1_X } / 

{[(Px+i)/(Px-or*+i} a) 

with 0 < X < 1 and P x = L was used to distribute the 
grid points along the X-direction. Three different grids 
(P x = 1.0, 1.01 or 1.05) were considered respectively in 
the computations. 

Within each duct cross section, the grid points were 
uniformly distributed along the 0-direction. Along the 
r-direction, the grid points were packed at the near duct 
wall locations (There were five grid points in Y + < 1 0 at 
the downstream location and the first grid point from the 
surface was at Y + = 0.79). Since we did not change the grid 
point distribution within each duct cross section, we only 
considered the effect of the X -directional grid packing on 
the computation. 

Initial Conditions . The inlet free stream properties 
were used as the initial conditions at each grid point to start 
the time iterations. The free stream had a total pressure 
14.7 psia, a total temperature of 520 °R, and a Mach 
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number of 0.34. These free stream conditions gave a bulk 
Reynolds number Re^ = 1.6xl0 6 . The Baldwin-Lomax 
eddy viscosity formula was used for turbulence modeling. 
The computation did not require turbulence initial 
conditions. 

Square Duct Flow 

A sketch of the flow configuration was shown in 
Fig. 3. We used the Cartesian coordinate system in the 
computation. The duct axis was in line with the inlet flow 
direction and the shape of the duct cross section did not 
change along the axial direction. Gessner 18 measured the 
flow properties within the near wall region. These existing 
experimental results could be used to verify the accuracy 
of the Navier-Stokes computations. The present 
computation was intended to predict some of the 
experimental results with bulk Reynolds number, Re^, = 
2.5x10 s . The existing experimental results indicated that 
the flow would be fully developed at X/D = 80. Thus, the 
computational domain was chosen to be a quadrant of the 
duct between X/D = 0 and X/D = 80. 

Boundary Conditions . Uniform velocity and uniform 
temperature were assumed at the duct inlet to simulate the 
flow with combined development of the hydraulic and 
thermal entry length. Near the duct inlet, the static pressure 
was assumed to vary linearly along the X-direction. Nonslip 
velocity conditions were the boundary conditions at the 
duct wall surface. The wall surface temperature was 
assumed to be uniform and the temperature was set at 1 . 1 
T m . The normal pressure gradients at the wall surface were 
set to zero. The flow symmetric conditions were the 
boundary conditions along the boundary at Y = D/2 and 
Z = D/2. Within the downstream duct cross-section (at 
X/D = 80), the X-direction gradients of the mean velocity 
components and the mean temperature were zero. They 
were the boundary conditions. However, a X-direction 
pressure gradient was imposed at the downstream location 
(X/D = 80). This pressure gradient was estimated according 
to the theoretical analysis 17 and the measurements of wall 
stresses. 18 

The Baldwin-Lomax eddy viscosity model and the 
Chien’s k-e turbulence model were studied respectively 
for the turbulence modeling in the computation. For the k 
ande computations, the gradients of k and e were set to 
zero at the upstream and the downstream boundaries. 
Noslip boundary conditions (k - 0 and e = 0) were used at 
the duct wall surface. Symmetry conditions were used at 
the planes of symmetry. 

Computational Grids . The grid dimensions were 
chosen to be 100 along the X-direction, 40 along the 
Y-direction, and 40 along the Z-direction. The grid 
generation technique in Ref. 10 was used to pack the grid 
points at the locations near the duct inlet and the duct wall 


surface. The grid point spacing along the X, Y, and Z 
directions were adjusted by specifying the packing 
parameters p x , P y and (5 Z . To study the grid effect on the 
computation, two different grids were used respectively to 
perform similar flow computations. 

One grid was generated with p x = 1.1, p Y = 
P z = 1.0075. The grid points were sparsely distributed at 
the near wall surface locations (There were 3 grid points 
in Y + < 10 at X/D = 80, Z/D = 0.5 and the first grid point 
from the surface was at Y + = 5.57). This grid was used in 
two similar computations with either the Baldwin-Lomax 
eddy viscosity formula or the Chien’ s k and e equations for 
turbulence modeling. Another grid was generated with 
P x = 1 . 1 , and P Y = P z = 1 .002 to increase the grid points 
in the near surface region (There were 5 grid points in 
Y + < 1 0 at X/D = 80, Z/D = 0.5 and the first grid point from 
the surface was at Y + = 1.84). This grid was used in the 
computation only with the Baldwin-Lomax eddy viscosity 
formula for turbulence modeling. 

Initial Conditions . We used two different turbulence 
modeling techniques in the square duct flow computations . 
Different initial conditions, corresponding to different 
turbulence modeling techniques, were used to start the 
time iterations. 

When the Baldwin-Lomax eddy viscosity formula 
was used in the computation, it used the upstream free 
stream properties as the initial conditions at each grid point 
location. The computed flow properties were the initial 
conditions for the square duct flow computation with the 
k-e turbulence model. The Proteus code could generate 10 
the initial k and e conditions. 

Circular to Rectangular Transition Duct Flow 

A sketch of the duct flow configuration was shown in 
Fig. 4. Details of the duct geometry and its setup in the 
wind tunnel experiment could be found in Ref 1 . We did 
not consider the actual shape of the wind tunnel inlet in the 
present computation. However, the circular section was 
extended upstream in the present computation. The 
extension was to account for the effect of the wind tunnel 
inlet shape on the flow development. For the present 
computation, the duct geometry consisted of an upstream 
circular section, a transition mid section, and a downstream 
rectangular section. The radius of the circular section was 
0.36 ft. The duct axis was in line with the inlet flow 
direction. The inlet flow had a Mach Number of 0.34, a 
total temperature of 520 °R and a total pressure of 1 4. 7 psia. 
These flow conditions gave a bulk Reynolds number. Re 
= 1 .6xl0 6 . We used the cylindrical coordinate system for 
the computation. 

Boundary Conditions . The computational domain 
was a quadrant of the duct between X/R = -6 and 
X/R = 7. The boundary conditions were specified at the 
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inlet, the duct surface, the symmetric plans, the duct center 
line, and the downstream boundary of the domain of 
computation. The boundary conditions were similar to 
what were previously described for the circular duct flow 
computation. A summary of the boundary conditions 
were described mathematically in the following: 

At the inlet location, X/R = -6, dp/dX = 0, 
U = U^, V = 0, W = 0, and T = T to . At the radial location, 
0 < r < R, along 0 = 0° direction, dp/d0 = 0, dU/d0 = 0, V 
= 0, dW/d0 = 0, and dT/d0 = 0. At the radial location, 0 < 
r < R, along = 90° direction, dp/d0 = 0, dU/d0 = 0, dV/d0 
= 0, W = 0, and dT/d0 = 0. At the location with r - 0 and 
0° < 0 < 90°, dp/dr = 0, dU/dr = 0, dv/dr = 0, dW/dr = 0, 
and dT/dr = 0. At the location with r = R and 0° < 0 < 90°, 
dp/dn = 0, U = 0, V = 0, W = 0, and T = Tw. At the 
downstream boundary location, X/R = 7, p/p^ = 0.985, 
dU/dX = 0, dV/dX = 0, dW/dX = 0, and dT/dX = 0. 

The downstream static pressure boundary condition 
(P^Poo = 0.985) was chosen from the experimental results. 1 
The purpose of the computation was to predict the surface 
heat transfer coefficients. Therefore, the surface tempera- 
ture of the experiment 6 was used as the surface temperature 
boundary condition at r = R. The experiment imposed a 
step increase in the wall surface temperature along the 
X-direction and the surface temperature was uniform 
along the 0-direction. Thus, the surface temperature 
boundary conditions for the present computation were 

T w = T 0 „. for X/D < -4, 

and 

T w = 1 . 1 97^, for X/D > -4. 

Computational Grids . The computations were 
performed with two different grids in order to study the 
grid effect on the computation. The Gridgen program 7 
was used to generate the grids. The grids had the same 
dimensions. The dimensions were 8 1 along the X-direction, 
50 along the r-direction, and 3 1 along the 0-direction. The 
upstream circular section and the transition mid section, 
respectively, had 3 1 *50*3 1 grid points. The downstream 
rectangular section had 21*50*3 1 grid points. Additional 
details of each grid were briefly described in the following 
sections. 

First Grid . The grid point spacing along the axial 
direction was uniformly distributed within each of the 
three different duct sections. Within a duct cross section, 
the grid points were packed at near duct surface locations 
(the grid point spacing was approximately the same as that 
of the circular duct flow computation). The grid point 
were uniformly distributed (A0 = constant) along the 


0-direction. Thus the grid orthogonality at the duct 
surface was maintained within the upstream circular section 
and along the side wall and the top wall center lines of the 
transition mid section and the downstream rectangular 
section. 

Second Grid . In the upstream circular section, the 
grid points fell on the duct cross sections. The grid points 
along the X-direction were packed at the locations near the 
duct entrance. Within each cross section, the grid points 
along the r-direction were packed at the locations near the 
duct axis and the duct wall surface (the near surface grid 
point spacing was approximately the same as that of grid 
1 ). At the same r location, the grid points were uniformly 
distributed in the 0-direction. Within the transition mid 
section, the grid orthogonality at the locations near the 
wall surface was optimized with Gridgen routines. The 
grid points were packed at the locations near the duct axis 
and the duct wall surface. Within the downstream 
rectangular section, the grid points fell on the duct cross 
sections. The grid points were equally spaced along the 
X-direction. The grid configuration within each cross 
section was the same as that at the downstream side 
(X/R = 3.75) of the transition mid section. 

Turbulence models . The Baldwin-Lomax eddy 
viscosity model was found effective (see the Results and 
Discussions section) for the computations of the internal 
flows within the circular duct and the square duct. This 
turbulence modeling technique was used for the circular to 
rectangular duct flow computation. 

Initial Conditions . The Proteus code was previously 
used to compute the turbulent developing circular duct 
flow properties. The computed circular duct flow properties 
were the initial conditions for the computation of the 
circular to rectangular duct flow. The Baldwin-Lomax 
eddy viscosity formula was used for turbulence modeling 
and the computation did not require turbulence initial 
conditions. 

Results and Discussions 

The previously described computational conditions 
were implemented to the Proteus computational code. The 
computations were proceeded to calculate the properties 
of the developing flows within the circular duct, the square 
duct and the circular to rectangular transition duct. The 
experience from the numerical computations indicated 
that the CFL value affected the starting of the computation. 
When the computation used the Baldwin-Lomax eddy 
viscosity, it could be started with a large CFL value 
(CF1 =10 was used here). A small CFL value (= 1) was 
required to proceed the computation when Chien’s k-E 
turbulence model was used in the computation. The 
computations converged within 45,000 time iterations. 
Some of the computed flow parameters, such as the duct 
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center-line velocity, the surface pressure, the skin friction 
factor and the surface heat transfer coefficient for each 
duct flow were described and compared with their values 
from theoretical analysis or experiments in the following 
sections. 

Circular Duct Flow 

For a developing cylindrical duct flow , the flow could 
be fully developed at downstream locations, X/R > 80. 17 
The present computational domain, L/R, along the 
X-direction was chosen to be either 1 00 or 140. The static 
pressure was also imposed as a boundary condition at the 
downstream location. The downstream pressure could be 
different at different L/R locations. The present 
computation was repeated with a few combinations of the 
domain configuration and the downstream pressure value. 
We investigated the effect of the domain and the 
downstream static pressure on the computation. 

A summary of the results of the X-directional 
variations of the center-line mean velocity, were 

shown in Fig. 5. The existing measurements of U^/U^ 
from experiment 17 were also plotted in Fig. 5 to compare 
with the present computed U c /U oo values. Different grid 
point packings (p x = 1.00, 1.01 or 1.05) along 
the X-direction did not affect the computation. A 

slight increase in the downstream pressure value (from 
0.945 to 0.960) did not change significantly the L^/U^ 
computations. With (i x = 1 .05, p d = 0.945p oc and different 
domain sizes (L/R = 100 or 140), the computations 
calculated approximately the same variations along 

the X-direction. All the computed L^/U^ values at near 
inlet locations (X/R < 30) agreed very well with the 
existing measurements. However, the computed 
values at X/R > 30 locations were found to be larger than 
their measured values. 

The computational results of the X-directional 
variations of the center-line static pressure, p^p^, were 
plotted in Fig. 6. The results showed that different static 
pressure level was computed when different 
computational domain or different downstream static 
pressure value was used in the computation. With L/R = 
140 or Pjj/p^ = 0.960 in the computation, it predicted an 
increase in the center-line static pressure within the entire 
duct flow region. The results showed that the static pressure 
varied linearly along the X-direction. 

The computational results of the skin friction factor, 
C f , were plotted in Fig. 7. The experimental results of the 
skin friction factors 17 were shown in this figure for 
comparison. Regardless of the difference in grid point 
packing, the computed skin friction factors at the near inlet 
locations (X/R < 30) agreed with the experimental values. 
The C f value decreased from a large value at the duct inlet 


to a small value at X/R = 30 location and then C f increased 
slightly at large X/R locations. When uniform grid point 
spacing (P x = 1 ) was used in the computation, the computed 
C f values oscillated at the downstream locations 
(X/R > 50). The computed values were well correlated 

in term of X/R. The skin friction factor of fully developed 
turbulent circular duct flow was correlated 17 in term of the 
bulk Reynolds number, 

Cf = [41o gl o{Re b /(4.531og IO Re b -3.82)}] _2 (2) 

The present duct flow had a bulk Reynolds number, 
Re b = 1.6xl0 6 . Equation (2) predicted C f = 2.70xl0 -3 for 
a fully developed flow. This skin friction factor value was 
also shown in Fig. 7 and compared with the present 
computational results of the downstream skin friction 
factors. The present computations predicted higher 
downstream skin friction factor value than it was computed 
with Eq. (2). 

The computed surface heat transfer coefficients, h, 
were plotted in Fig. 8. The heat transfer coefficient was 
large at the duct inlet location and the h value decreased 
gradually at the downstream locations. The size of the 
domain of computation and the downstream static pressure 
value seemed to have little effect on the computation of the 
surface heat transfer coefficient. However, the computed 
h values at the downstream locations oscillated when the 
grid points were distributed uniformly (f$ x = 1 ) along the 
X-direction. The h values were well correlated in term of 
X/R values. It was established 16 that the Nusselt number, 
Nu (= HD/K), and the bulk Reynolds number of a fully 
developed turbulent circular duct flow were correlated 
with 

Nu = 0.021 Pr° 6 Re b 8 , for 10 4 < Re b < 10 6 (3) 

Substituting the present bulk Reynolds number, 
Re b = 1 .6x1 0 6 , in Eq. (3), it predicted a Nusselt number of 
1.5xl0 3 . This Nusselt number corresponded to a heat 
transfer coefficient, h = 755. This value was also shown in 
Fig. 8 and it was larger than the computed h values at the 
downstream locations. The present computation indicated 
that, with combined hydraulic and thermal entry length, 
the circular duct flow had large thermal entry length. This 
is similar to the analysis of the thermal entry length with 
constant heat transfer rate. 16 

Square Duct Flow 

Existing experiment 18 indicated that a developing 
turbulent square duct flow was fully developed at 


6 



X/D =80. For a fully developed square duct flow, 
analysis 17 showed that the averaged wall shear stress, f w , 
and the X-di recti on pressure gradient were related with 

dp/dx = -4 t w /D (4) 

Based on the experimental results of the averaged 
shear stress, 1 8 we used Eq. (4) to calculate the X-direction 
pressure gradient. The pressure gradient had a value of 
-0.045 lb/ft 2 . This pressure gradient value was used as the 
downstream pressure boundary condition at X/D = 80 for 
the present computation. Some computational results of 
the duct center-line velocity, the skin friction factor and 
the surface heat transfer coefficient are described in the 
following. 

The X-directional variations of the duct center-line 
velocity, U C /U M , were plotted in Fig. 9. The computational 
results due to different grids and different turbulence 
models were compared with the existing experimental 
results. 19 All the results indicated that U c /U oe increased to 
a large value at X/D = 40 and then U c /U oo decreased slightly 
to a constant value at the downstream locations. The 
computational results showed a grid dependence when the 
Baldwin-Lomax eddy viscosity formula was used for 
turbulence modeling. The computation with p x = 1 .1 and 
P y = P z = 1.0075 predicted better the measured l^/U^ 
valued than the computation with (3 X = 1.1 and p Y = 
P z = 1.002. With the Chien’s turbulence model and 
p x = 1 . 1 and P Y = P 2 = 1 .0075 grid in the computation, it 
could predict the measured U c /U eo at most X/D locations. 
The computation did not predict better the measured U c / 

values than the computation with the Baldwin-Lomax 
eddy viscosity formula for turbulence modeling. Thus, the 
Baldwin-Lomax turbulence model was effective for the 
present computation to calculate the center-line velocity. 

The computational results of the skin friction factor 
variations along the perimeter of the duct cross section at 
X/D = 80 were plotted in Fig. 1 0. The results with different 
grid configurations and different turbulence models were 
compared with existing experimental results. 18 These 
results indicated that C f had a very small value at the 
comer location and the C f value increased rapidly at the 
near comer locations (Z/D < 0.1) to a large value at 
Z/D = 0.25 location. The C f value then remained 
approximately the same at the other peripheral locations. 
With the Baldwin-Lomax eddy viscosity formula and P x 
= 1.1 and p y = p z = 1 .002 grid, the computed C f values 
agreed very well with the C f measurements at all Z/D 
locations. With the Baldwin-Lomax eddy viscosity formula 
and p x = 1 . 1 and P Y = P z = 1 0075 grid in the computation, 
it calculated small friction factors at the near comer 
locations. With the Chien ’ s turbulence model and p x = 1. 1 
and P Y = P z = 1 0075 grid in the computation, it computed 


high skin friction factors at all peripheral locations. With 
grid point packing at the near wall locations and the 
Baldwin-Lomax eddy viscosity formula for turbulence 
modeling, the computation could predict the experi mental 
results of the skin friction factors. 

The computational results of the surface heat transfer 
coefficients along the perimeter of the duct cross section 
at X/D = 80 were plotted in Fig. 1 1 . The computational 
results due to different turbulence models or different 
computational grids were compared in this figure. The 
results indicated that very small heat transfer coefficients 
were calculated at comer locations. With the Baldwin- 
Lomax eddy viscosity formula and either one of the two 
different grids in the computation, it calculated 
approximately the same level of the surface heat transfer 
coefficients. With the Chien’s k-e turbulence model and 
p x =1.1, and (J Y = p z = 1 .0075 grid in the computation, it 
calculated high surface heat transfer coefficients. 

Circular to Rectangular Transition Duct 

As it was previously described in the Computation 
section, the duct inlet was circular and the circular duct 
section was extended upstream to account for the actual 
wind tunnel inlet in the experiment. 6 A downstream static 
pressure, pyp^ = 0.985, was imposed at the downstream 
boundary of the domain of computation. The computation 
calculated a uniform static pressure (p^p^ = L00) at the 
duct inlet. The computed Mach number profile along the 
direction normal to the duct surface at X/R = -2 was 
plotted in Fig. 12. This Mach number profile was also 
compared with its experimental results. 1 The comparison 
showed good agreement between the results from the 
computation and the experiment. Some computational 
results of the surface pressure, the skin friction factor, and 
the heat transfer coefficient were described in the following 
sections. 

The computational results of the static pressure 
variations along the top wall center line and the side wall 
center line were plotted in Fig. 13. The computational grid 
induced very small difference between the static pressures 
along the wall center lines. The difference is less than 
0.01 p^. Neglecting the grid effect on the static pressure 
computations, we could observe general trend of the static 
pressure variations along the wall center lines. Along the 
side wall center line locations, the pressure decreased and 
then increased within most part of the transition section. 
The pressure decreased rapidly within the downstream 
side of the transition section and the pressure maintained 
the same value along the side wall center line of the 
rectangular section. Along the top wall center line, the 
pressure increased and then decreased within most part of 
the transition section. The pressure increased rapidly at 


7 



the end of the transition section and the pressure maintained 
the same value along the top wall center line of the 
rectangular section. The pressure levels were the same 
along the top wall center line and the side wall center line 
of the downstream rectangular section. Experiment 1 
measured the pressure along the top wall center line. The 
measured values were plotted in Fig. 13 to verify the 
accuracy of the pressure computation. The 
computed pressure values were found to agree very well 
with their experimental values. 

The computed skin friction factor, C f , at locations 
along the top wall and the side wall center lines were 
shown in Fig. 14. The computational grids only induced 
difference in the predictions of the center-line C f values 
within the upstream circular duct section. However, the C f 
variation within the circular section was similar to the C f 
variation of the developing circular duct flow (Fig. 7). The 
top wall center-line C f values were the same as the side 
wall center-line C f values. This trend was reversed within 
the upstream region of the transition section. In that 
region, the top wall center-line C f values were lower than 
the side wall center-line C f values. At the downstream 
locations, the top wall center-line C f were higher than the 
side wall center- line C f . Along the side wall center-line, C f 
decreased to a small value and then increased within the 
transition section. The C f values were approximately the 
same within the rectangular section. Along the top wall 
center-line, C f decreased at the upstream locations of the 
transition section. Then C f increased to a large value 
within the transition duct section. Again C f decreased 
rapidly at locations near the end of the transition section. 
C f maintained approximately the same values within 
rectangular section. 

The surface heat transfer coefficient, 

H = Qw /(T w - T.) = -KdT/dn| n=0 /(T w - T.) 

was computed at the locations along the top wall center line 
and the side wall center line. The computed surface heat 
transfer coefficients with the first grid in the computation 
were plotted in Fig. 15. In the computation, a step increase 
in the surface temperature boundary condition was imposed 
at X/R = -4 location (Fig. 4). At this location, H increased 
rapidly to a large value. H then decreased within the circular 
section and the H value at the side wall center line was about 
the same as the H values along the top wall center line. Along 
the side wall center line of the transition section, H decreased 
to asmall value and then increased rapidly to aconstant value 
within the rectangular section. Along the top wall center line, 
H increased at the upstream locations of the transition 
section. The H value then decreased slightly and maintained 
at a constant level at the downstream locations. In the 
rectangular section, the top wall center line had a larger H 


value than that of the side wall center line. Figure 1 6 showed 
the center line H values which were calculated with the 2nd 
grid in the computation. The computed H along the center 
lines were similar to the results (Fig. 1 6) when the 1st grid 
was used in the computation. 

The experimental results 6 of the center-line heat 
transfer coefficients were plotted in Figs. 15 and 16 to 
verify the accuracy of the present heat transfer computation. 
The heat transfer coefficient of a fully developed turbulent 
circular duct flow 16 was also indicated in these figures. 
The computation predicted very well the experimental 
heat transfer coefficients at the end of the circular section. 
The heat transfer coefficients at X/R = 0 were very close 
to the H value of a fully developed turbulent circular duct 
flow. The present computation could predict the variation 
of the measured H value along the side wall center line. In 
the transition section, the computed top wall center line 
heat transfer coefficients were lower than the experimental 
H values. Using the 2nd grid in the computation, it 
improved the prediction of the top wall center-line heat 
transfer coefficients. 

The computed H values at the locations along the 
perimeters of several duct cross sections were presented in 
a polar plot format in Fig. 17. Since the computation with 
the 2nd grid improved the prediction of the H values along 
the center line locations, the peripheral H variations due to 
2nd grid in the computation were shown in Fig. 1 7. Some 
experimental H values were also plotted in this figure for 
comparisons with the computed H values. At X/R = - 1 , the 
duct had a circular cross section and the computation 
calculated approximately the same H values at all peripheral 
locations. The computed H values agreed well with their 
experimental values. At X/R = 0, the duct had a circular 
cross section and the computation predicted a decrease in 
the H values at large 0 locations. The computational result 
of the H variation was similar to the result from the 
experiment. Within the transition section (at 
X/R = 1 .0, 1 .5 and 2), the computed H values decreased on 
the side wall surface and the H increased rapidly at the 
comer locations. The top wall surface H values were 
higher than the side wall surface H values. The computation 
could predict the experimental results of the peripheral 
heat transfer coefficient variation. The computed H values 
at all peripheral locations were large at the end of the 
transition duct section (X/R = 3.5) and the computation 
predicted accurately the experimental heat transfer 
coefficients at the top wall surface. The results in Fig. 17 
also indicated that the transition section had large heat 
transfer coefficients at the duct comer locations. This 
observation was contrary to the results of the square duct 
computation which predicted very small heat transfer 
coefficients at the square duct comer locations. This could 
be due to the generation of the vortices at the transition 
duct comer locations. 1 
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Summary 

The objective of this research was to study the 
computation of the surface heat transfer coefficients of a 
transition duct flow. The transition duct flow had combined 
development of hydraulic and thermal entrance length. 
The axis of the transition duct was in line with the inlet 
flow direction. The inlet flow was subsonic (Mach number 
= 0.34). The bulk Reynolds number was 1 .6x 1 0 6 . The duct 
surface temperature was uniform. The duct cross section 
changed gradually from a circular shape at the upstream 
locations to a rectangular shape at the downstream 
locations. 

The existing Proteus three-dimensional Navier- Stokes 
numerical computational code was used for the 
computations of the heat transfer coefficients. The Gridgen 
was used to generate the computational grids. Three major 
computational procedures were performed in this heat 
transfer study. The Navier-Stokes code was first used to 
compute the properties of the developing flows within a 
circular duct and a square duct. These computations served 
as the base-line verification of the effect of the 
computational grid, the Baldwin-Lomax eddy viscosity 
formula and the Chein’s k-e turbulence model on the 
computation. The code was then used to compute the 
surface heat transfer coefficients of the transition duct 
flow. 

The present computation, with the Baldwin-Lomax 
eddy viscosity formula for turbulence modeling, predicted 
approximately the experimental results of the variations 
of the center-line velocity and the skin friction factor 
within the entrance region of a circular duct flow. The 
computation required grid point packing in the duct 
entrance in order to predict simultaneously the existing 
results of the center-line velocity, the skin friction factor 
and the surface heat transfer coefficient. The grid packing 
along the axial-direction had little effect on the prediction 
of the center-line velocity. 

The square duct flow computation, with proper grid 
point packing and the Baldwin-Lomax eddy viscosity 
formula for turbulence modeling, predicted the 
measurements of the center-line velocity within the 
entrance region. With a few grid points in the near wall 
region, the computation also predicted the measurements 
of the skin friction factors along the perimeter of the 
downstream duct cross section. The computed surface 
heat transfer coefficients had very small values at the near 
comer locations. With the Chien’s k-e turbulence model 
in the computation, it improved the predictions of the 
center-line velocities at the downstream location and 
calculated large skin friction factors and the surface heat 
transfer coefficients at downstream locations. 

The computations of the circular duct flow and the 
square duct flow indicated that the Baldwin-Lomax eddy 


viscosity formula was effective for turbulence modeling 
in the heat transfer computation. We only used this 
turbulence model in the transition duct computation. The 
computations were performed with two different grid 
configurations. The computation predicted the patterns 
of experimental surface heat transfer coefficient variations. 
The computed heat transfer coefficients at the upstream 
locations agreed very well with their experimental values 
and the computed heat transfer coefficients at the 
downstream locations were lower than their 
measurements. Different grid configuration induced small 
differences in the computational results of the surface 
pressure, the skin friction factor and the surface heat 
transfer coefficient. The computed surface heat transfer 
coefficients had large values at the locations near the 
comer of the downstream rectangular duct section. 

The results in this report were obtained with the 
straight forward application of the existing Proteus three- 
dimensional Navier-Stokes computational code. The 
present computation used many default options of the 
numerical computational methods in the code. Some 
numerical schemes, such as the thin-layer option, various 
artificial viscosity models and different convergence 
criteria were not studied for their effects on the 
computations. The Proteus code provided great flexibility 
to specify the boundary conditions for the Navier-Stokes 
computations of an internal flow properties. We took the 
flexibility to specify the boundary conditions for the 
present computations. Theoretical analysis of the boundary 
conditions plays an important role in the numerical 
computations of subsonic flow properties. A detailed 
study of the transition duct flow boundary conditions is 
required to pursue further the flow computation. The 
present research considered only simple internal duct 
flows. The boundary conditions were derived from simple 
analyses and experimental results of the flow properties. 
The comparisons between the computed heat transfer 
coefficients and their measurements indicated that the 
Navier-Stokes computation could predict the heat transfer 
coefficients of a transition duct flow. The experience 
learned from the research was that a Navier-Stokes 
computation of the surface heat transfer rates of a transition 
duct, including the influence of boundary conditions, grid 
configuration and turbulence modeling, could be a 
formidable task. 
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